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Abstract. The equilibrium distributions of a Markovian model describ- 
ing the interaction of several classes of permanent connections in a net- 
work are analyzed. It has been introduced by Graham and Robert [5]. 
For this model each of the connections has a self-adaptive behavior in 
that its transmission rate along its route depends on the level of conges- 
tion of the nodes on its route. It has been shown in that the invariant 
distributions are determined by the solutions of a fixed point equation in 
a finite dimensional space. In this paper, several examples of these fixed 
point equations are studied. The topologies investigated are rings, trees 
and a linear network, with various sets of routes through the nodes. 



1 Introduction 

Data transmission in the Internet network can be described as a self-adaptive 
system to the different congestion events that regularly occur at its numerous 
nodes. A connection, a TCP flow, in this network adapts its throughput accord- 
ing to the congestion it encounters on its path: Packets are sent as long as no 
loss is detected and throughput grows linearly during that time. On the contrary 
when a loss occurs, the throughput is sharply reduced by a multiplicative fac- 
tor. This scheme is known as an Additive Increase and Multiplicative Decrease 
algorithm (AIMD). 

Globally, the TCP protocol can be seen as a bandwidth allocation algorithm 
on the Internet. From a mathematical modelling perspective, the description 
is somewhat more difficult. While the representation of the evolution of the 
throughput of a single TCP flow has been the object of various rigorous works, 
there are few rigorous studies for modelling the evolution of a large set of TCP 
connections in a quite large network. 
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A possible mathematical formulation which has been used is via an optimiza- 
tion problem: given K classes of connections, when there are Xk connections of 
class k € {!,..., K}, their total throughput achieved is given by Afe so that the 
vector (Afc) is a solution of the following optimization problem 



where A is the set of admissible throughputs which takes into account the capac- 
ity constraints of the network. The functions (C/fc) are defined as utility functions, 
and various expressions have been proposed for them. See Kelly et al. 7 , Mas- 
soulie [5] and Massoulie and Roberts [TUJ. With this representation, the TCP 
protocol is seen as an adaptive algorithm maximizing some criterion at the level 
of the network. 

A different point of view has been proposed in Graham and Robert It 
starts on the local dynamics of the AIMD algorithm used by TCP and, through 
a scaling procedure, the global behavior of the network can then be described 
rigorously. It is assumed that there are K classes of permanent connections 
going through different nodes and with different characteristics. The loss rate of 
a connection using a given node j, 1 < J < J, is described as a function of the 
congestion Uj at this node. The quantity Uj is defined as the (possibly weighted) 
sum of the throughputs of all the connections that use node j. The interaction 
of the connections in the network is therefore expressed via the loss rate at each 
node. 

It has been shown in Graham and Robert [5] that under a mean-field scal- 
ing, the evolution of a class k connection, I < k < K, can be asymptotically 
described as the unique solution of an unusual stochastic differential equation. 
Furthermore, it has also been proved that the equilibrium distribution of the 
throughputs of the different classes of connections is in a one to one correspon- 
dence with the solution of a fixed point equation (£) of dimension J (the number 
of nodes). 

Under "reasonable" conditions, there should be only one solution of {£) and 
consequently a unique stable equilibrium of the network. Otherwise this would 
imply that the state of the network could oscillate between several stable states. 
Although this is mentioned here and there in the literature, this has not been 
firmly established in the context of an IP network. It has been shown that multi- 
stability may occur in loss networks, see Gibbens et al. [4] and Marbukh [8] or 
in the context of a wireless network with admission control, see Antunes et 
al. [1]. Raghunathan and Kumar [12] presents experiments that suggest that a 
phenomenon of bi-stability may occur in a context similar to the one considered 
in this paper but for wireless networks. 

It turns out that it is not easy to check in practice whether the fixed point 
equation {£) has a unique solution or not. The purpose of this paper is to in- 
vestigate in detail this question for several topologies. The paper is organized as 
follows. Section [2] reviews the main definitions and results used in the paper. In 
addition, a simple criterion for the existence of a fixed-point solution is given. 
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Section [3] presents a uniqueness result for a tree topofogy under the assump- 
tions that ah connections use the root. Section |4] considers a linear network. 
Section [5] studies several scenarios for ring topologies and a uniqueness result 
is proved for connections going through one, two, or all the nodes. Two main 
approaches are used to prove uniqueness: monotonicity properties of the network 
and contraction arguments. 

A general conjecture that we make is that when the loss rates are increasing 
in the level of congestion, this should be sufficient to imply the uniqueness of the 
equilibrium in a general network (together with regularity properties perhaps). 

2 A Stochastic Fluid Picture 

In this section, a somewhat simplified version of the stochastic model of inter- 
acting TCP flows of Graham and Robert [5] is presented. 

The case of a single connection 

Ott et al. [TT] presents a fluid model of a single connection. Via scalings with 
respect to the loss rate, Dumas et al. [3] proves various limit theorems for the 
resulting processes. The limiting picture of Dumas et al. 3J for the evolution of 
the throughput of single long connection is as follows. 

If the instantaneous throughput at time t of the connection is W{t), this 
process has the Markov property and its inflnitesimal generator is given by 

n{f){x) = af{x) + (3x{f{rx) - }{x)) (1) 

for / a C^-function from M+ to M. For t > 0, the quantity W{t) should be 
thought as the instantaneous throughput of the connection at time t. 

The Markov process {W{ty) increases linearly at rate a. The constant a is 
related to the distance between the source and the destination. It increases pro- 
portionally to the round trip time RTT, typically 

" Ci+RTT' 

for some constants Co and Ci. 

Given W{t) = x, the process {W{t)) jumps from x to rx [r is usually 1/2) at 
rate (3x. The expression I3x represents the loss rate of the connection. Of course, 
the quantities a, /3 and r depend on the parameters of the connection. 

The density of the invariant distribution of this Markov process is given in 
the following proposition. It has been analyzed in Ott et al. [TT] at the fluid level 
and by Dumas et al. [3] , see also Guillemin et al. (6] . The transient behavior has 
been investigated in Chafai et al. [2] . 

Proposition 1. The function 

H (w)- V — e-'"-""'"'/2 w>0 (2) 
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with p — a/ [3, is the density of the invariant distribution of the Markov process 
(W{t)) whose infinitesimal generator is given by Equation Furthermore, its 
expected value is given by 

/ wHr,p{w) dw = J—\{ _ . (3) 

•'0 n=l 

A Representation of Interacting Connections in a Network 

The network has J > 1 nodes and accommodates K > 1 classes of permanent 
connections. For 1 < k < K, the number of class k connections is Nk > 1, and 
one sets 

N ={Nu...,Nk), and \N\= Ni + --- + Nk. 

An allocation matrix A = [Ajk, 1<J<J, l<fc< K) with positive coefficients 
describes the use of nodes by the connections. In particular the route of a class k 
connection goes through node j only if Ajk > 0. In practice, the class of a 
connection is determined by the sequence set of nodes it is using. 

If Wn^k > is the throughput of the nth class k connection, I < n < Nk-, 
the quantity AjkWn,k is the weighted throughput at node j of this connection. 
A simple example would be to take Ajk = 1 or depending on whether a class k 
connection uses node j or not. The total weighted throughput Uj of node j by 
the various connections is given by 

""j = ^ ^ AjkWn,k- 
k=l n=l 

The quantity uj represents the level of utilization/congestion of node j. In par- 
ticular, the loss rate of a connection going through node j will depend on this 
variable. 

For I < k < K, the corresponding parameters a and (3 of Equation ^ for 
a class k connection are given by a non-negative number Uk and a function Pk '■ 
M."^ R+, so that when the resource vector of the network is u = {uj, I < j < J) 
and if the state of a class k connection is vuk- 

— Its state increases linearly at rate Ofc. For example Ok ~ ^/Rk where Rk 
is the round trip time between the source and the destination of a class k 
connection. 

— A loss for this connection occurs at rate WkPk{u) and in this case its state 
jumps from Wk to rkWk- The function Pk depends only on the utilization of 
all nodes used by class k connections. In particular, if a class k connection 
goes through the nodes , j2 , ■ • ■ , i one has 

Pkiu) = PkiUj,,Uj^,...,Uj,^). 

A more specific (and natural) choice for Pk would be 



Pk{u) = <5fc + fj, {ujj + ipj^ (mjJ + • • • + ifjj^ (wj,^ ), (4) 
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where Lpj^ {x) is the loss rate at node j£ when its congestion level is a; > 0, 
and 5k is the loss rate in a non-congested network. Another example is when 
the loss rate (3k depends only on the sum of the utilizations of the nodes 
used by class k, i.e., 

Pk{u)=l3k{j^uj}j. (5) 



Asymptotic behaviour of typical connections 

If {Wn,k{t)) denotes the throughput of the nth class k connection, 1 < n < iVfc, 
then the vector 

{W{t)) = {[{WnAi)), l<k<K,l<n<NK],t>Q) 

has the Markov property. As it stands, this Markov process is quite difficult to 
analyze. For this reason, a mean field scaling is used to get a more quantitative 
representation of the interaction of the flows. More specifically, it is assumed 
that the total number of connections ||A^|| goes to infinity and that the total 
number of class k connections is of the order pfc||A^||, where pi + ■ ■ ■ + pk = 1- 

For each 1 < k < K, one takes a class k connection at random, let Uk be 
its index, I < Uk < Nk- The process {Wn^,k{t)) represents the throughput of a 
"typical" class k connection. It is shown in Graham and Robert that, as ||iV|| 
goes to infinity and under mild assumptions, the process [{Wn^^kit)), 1 < fc < if] 
converges in distribution to {W{t)) = [{Wk{t)), ^ < k < K], where the processes 
{Wk{t)), for 1 < fc < if, are independent and, for 1 < /c < if, the process 
{Wk{t)) is the solution of the following stochastic differential equation, 

dWkit)^akdt-{l-rk)Wkit-) J l[o<z<w,(t-)04u^it))}^kidz,dt), (6) 

with Uyy{t) = 1 < j < J) and, for 1 <j < J, 

K 

fe=i 

where (A/fc, 1 < A; < K) are i.i.d. Poisson point processes on with Lebesgue 
characteristic measure. 

Because of the role of the deterministic function {u-^^j^) in these equations, 
the Markov property holds for this process but it is not time-homogeneous. The 
analogue of the infinitesimal generator f2k^t is given by 

72fc^t(/)(x) = akf'ix) + xl3k{u^{mf{rkx) - f{rk)). 

The homogeneity holds when the function {u-^{t)) is equal to a constant u* , 
which will be the case at equilibrium. In this case a class k connection behaves 
like a single isolated connection with parameters a = ak and /? = j3k{u*). 
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The Fixed Point Equations 

The following theorem gives a characterization of the invariant distributions for 
the process {W{t)). 

Theorem 1. The invariant distributions for solutions {W{t)) of Equation ([6]) 
are in one-to-one correspondence with the solutions u G M.'j_ of the fixed point 
equation 

K 

uj = AjMu), 1 < J < J, (7) 

A,-l 

where 

// u* is such a solution, the corresponding invariant distribution has the density 
w — > Y[k=i -^rk.pki'^^k) on R-f , where pk — Ok/Pkiu*) and H^.p is defined in 
Proposition [7J 

The above theorem shows that if the fixed point equation ([7]) has several solu- 
tions, then the limiting process {W{t)) has several invariant distributions. Simi- 
larly, if equation ^ has no solution then, in particular, {W{t)) cannot converge 
to an equilibrium. These possibilities have been suggested in the Internet liter- 
ature through simulations, like the cyclic behavior of some nodes in the case of 
congestion. 

Under mild and natural assumptions, such as the loss rate being non- 
decreasing with respect to the utilization in the nodes, we show that for some 
specific topologies there exists a unique fixed point. We believe that such a 
uniqueness result will hold, in fact, for any network in general (under suitable 
regularity properties on the functions fc = 1, . . . , K). Before proceeding to the 
examples, we first present an existence result that holds for a general network. 

2.1 An Existence Result 

In this section the existence of a solution to the fixed-point equation ([7]) is proved 
for a quite general framework. 

If u is a solution of Equation ([7]) and Zfe = 4>k{u), I < k < K , then the vector 
z = (zk) satisfies the relation u — Az, i.e., Uj — Ajizi + Aj2Z2 + • ■ • + Ajkzk-, 
1 < j < J, and as well 

z = 'P{z)'^= {(t)k{M,'^<k<K). (9) 
The proposition below gives a simple criterion for the existence of a fixed point. 

Proposition 2. // the functions u j3k{u), 1 < k < K, are continuous and 
non- decreasing, and if there exists a vector z'^'^^ G R;^ such that the relations 

< oo. 
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hold coordinate by coordinate, then there exists at least one solution for the fixed 
point Equation ([9]) and therefore also for Equation ([7| • 

Proof. Define the sequence z'"^ = ^(z^""^'), n = 1,2, . . .. From z^"^ < ^(z(°)) 
and 2*°) < ^(<?>(2(°))), it follows that < z'^^ and z^") < z^^K Since the 
function ^ is non-increasing, one gets that the relation 

< < . . . < z(2") < . . . < < . . . < < ^(1) 

holds. Hence, there are z*, z* G M;^ such that 

lim z(2«) ^ and lim z'^"-!) = 

with z, < z*. Since z'^") = ^(z^^"-!)) and z(2"+i) = <?(z(2")), by continuity we 
also have that z* = <P{z^) and z* = <?(z*). 

Define the set = {z : z, < z < z*}. Note that z* < z^^) = ${z^°^) < oo, 
hence is bounded. In addition, for z Cz D, 

z* =${z*)<^{z) <^{z,)^z*, 

since the function ^ is non-increasing. One can therefore apply Brouwer fixed 
point theorem to ^ restricted to the compact convex set D, and conclude that 
D contains at least one fixed point of the function ^. The proposition is proved. 

The conditions of Proposition [2] trivially hold when f3k is non-decreasing and 
/3fe(0) > for all fc = 1, ... .if, since then ${0) < oo, < <?(0) and < ^^{0))- 
In particular, when the function Pk is given by ([3]) , 5^ > is a sufficient condition 
for the existence of a fixed point. 

3 Tree topologies 

We consider a finite tree network. A connection starts in the root and then follows 
the tree structure until it leaves the network at some node. The set of routes is 
therefore indexed by the set of nodes, i.e., a connection following route G € T 
starts in the root and leaves the tree in node G. 

The tree can be classically represented as a subset T of U„>oN" with the 
constraint that if G = (gi,...,gp) € T, then, for 1 < £ < p, the element 
H = {gi, . . . , g^) is a node of the tree as well. In addition, node H is the g^th 
child of generation (level) £ and the ancestor of G for this generation. One writes 
H C G in this situation and H \- G when £ — p — 1, i.e., when G is a daughter 
of H. The quantity U[h.g] denotes the vector {up, P : H C P C G). The root of 
the tree is denoted by 0. Assume that Ahg = 1 if route G uses node H, and 
otherwise. Equation ([7]) writes in this case, 

UH = X! '?^g("[0.G]), H eT, 

G&T,HCG 
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Fig. 1. Tree with connections starting at root node 

which is equivalent to the recursive equations 

UH = <I>h{u[i!im]) + X! ■"<^' HeT. (10) 

G£T,H\-G 

Proposition 3. If the functions [3h, H Cz T, are continuous and non- decreasing, 
then there exists a unique solution for the fixed point equation 

Proof. Let H he a, maximal element on T for the relation C, i.e., H is a leaf, 
and denote by P{H) the parent of node H. Equation (fTO|) then writes 

UH = (j)H{u[(ii^p(H)],UH)- (11) 

The function (pn being non-increasing and continuous, for a fixed vector U[0,p(_f/)] , 
there exists a unique solution uh = (■W[0_p(jy)]) > to the above equation. 
Furthermore, the function W[0^p(i/)] — > Fh{u[(Iip(^h)]) is continuous and non- 
increasing. For such an H, for H' — P{H), Relation (flOl) can then be written 
as 

UH' = 't'H'{u[$^P{H')],UH') + ^ Fg{u[$^p(^h')],Uh')- 

GeT,H'hG 

Since (j)H' and Fg, with G a leaf, are non- increasing and continuous, there exists 
a unique solution uh' — Fh'{u[id,p(h')]) ^ and the function 

"[0,P(H')] Fh'{u[II,^P{H')]), 

is continuous and non-increasing. By induction (by decreasing level of nodes), 
one obtains that a family of continuous, non-increasing functions Fq, G G T, 
G ^ 0, exists, such that, for a fixed vector U[0 p(g)], uq — ^g(w[0.p(g)]) is the 
unique solution of 

ug = 0g("[0,p(g)],'"g) + X! Pg'{uiid^p{g)],ug)- 

G'er,GhG' 

Equation PH)) at the root then writes 

m ^ him) + ^ Faiudi), 



9 



and this equation has a unique solution U0. Now, one defines recursively (by 
increasing level of nodes) 

UG = Fg{u[{/>,p(g)]), GeT. 
Then clearly {ug, G E T) satisfies Relation ([TOl) and is the unique solution. 

4 Linear topologies 

In this section we consider a linear network with J nodes and K ~ J + 1 classes 
of connections. Class j connections, I < j < J, use node j only, while class 
connections use all J nodes. Assume Ajk = 1 if class k uses node j, and 
otherwise. Equation ([7]) is in this case 

Uj = 00 (w) + (bj{uj), 1 < J < J, (12) 

with u = (ui, . . . , uj). 



1 


H 




>- -4 












5- 



Fig. 2. A linear network with J nodes and K = J + 1 classes of connections 



Proposition 4. // the functions, /3fc,0 < k < J , are continuous and non- 
decreasing, then there exists a unique solution for the fixed point equation ([7]). 

Proof. Let 4>j{x) — x — (j)j{x), x S M, which is continuous and non-decreasing. 
Hence, can be rewritten as 

u, = {Uu)) = ( ^= ) ,^<J<J- (13) 
\^(3o{u)) 

for some constant a^, see Equation ([5]). In addition, define the function ipj{x) = 
(f'J^ {oiq/ \/x), a: G M, which is continuous and non-increasing. ^^From (|13p we 
obtain the relation 

/3o(w) = /?o(^i(/3o(u)),...,Vj(/3o("))). 

Since (3q is non-decreasing and ij^j is non-increasing, the fixed point equation 
(3 = /?o(V'i(/3), ■ • ■ , has a unique solution (3* > 0. Hence, the Relation ifTS]) 

has a unique fixed point, which is given by u* = (f)J^ (q;o/a//3*). 
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5 Ring topologies 

In this section, the topology of the network is based on a ring. Several situations 
are considered for the paths of the connections. 

Routes with two consecutive nodes 

It is assumed that there are J nodes and K — J classes of connections and class 
j e {1, . . . , J} uses two nodes: node j and j + 1. Assume Ajk = 1 if class k uses 
node j, and otherwise. Equation ([7]) is in this case 

Uj = 4>j-i{uj^i,u.i) + (f>j{uj,u.i+i), j = l,...,J. (14) 

For Uj = (j)j{uj^ "^+1)5 the above equation can be rewritten as follows 

Vj = (t^jiVj-i + y^^V] + j = 1, ■ ■ • , ^- (15) 




Fig. 3. Routes with two consecutive nodes 



Proposition 5. If the functions [3k, 1 < A: < K, are continuous, non- decreasing 
and satisfy the assumptions of Proposition then there exists a unique solution 
for the fixed point equation ([7]) . 

Proof. ^From Proposition [5] we have that Equation (fT5|) has at least one fixed 
point solution. Let x = (xj : j = 1, 2, . . . , J) and y — {yj : j — 1,2, ... , J) both 
be fixed points. 

If the relation yj < Xj holds for all j = 1, . . . , J, then the inequality 

(t^jiVj +yj-i,y]+i +yj) = % < = ^^(xj +a;j-i,a;j+i +2;^), 

and the fact that the function (pj is non-increasing, give directly a contradiction. 

Consequently, possibly up to an exchange of x and y, one can assume that 
there exists m € {1,...,J} such that pm < Xm and j/m+i > Xm+i- Define 
Cj = Xj — yj and dj — yj — Xj. Hence, Cm > and dm+i > 0. Without loss of 
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generality, it can be assumed that the classes are ordered such that dm-i < dm+i- 
Since the function 0™ is non-increasing, and 

we have that cither 

and/or j/,„ + y„i+i > x 

i.e., dm^i > Cm and/or d^+i ^ Cm- Because dm-i < c^m+i, then, necessarily, 
dm^i > Cm > 0. Hence 

> Xm+l — 4>m+l{Xm+l + a^m, a:^m+l + 

Since 4>m+i is non-increasing, one has ym+i + ym-1-2 < Xm+i + a;,„+2 and conse- 
quently d„j+i < C„i+2- 

^From < Cm < dm+i < Cm+2, we obtain Xm+2 > ym+2, which, using the 
same steps as before, implies Cm-t-2 < dm+i- In particular, by induction it can be 
concluded that 

Cj < dj+i < Cj+2 < dj+3, for all j = 1, . . . , J, 

where the indices j + 1, 7 + 2, and j + 3 are considered as modulo J. This implies 
that Cj — dj — c, for all j = 1, . . . , L, and hence yj + j/j-i = Xj + i.e., 

Vj = + Vj-i^Vj + Vj+i) = -^jCa^j + Xj-i,Xj + .Tj+i) = Xj, 

for i — 1, . . . , J . We can conclude that the fixed point is unique. 

The rest of this part will be devoted to a contraction argument that can be 
used to get a unique solution to the fixed point equation. 

Proposition 6. // the functions (3k, 1 < fc < K, are Lipschitz, continuous 
dijferentiable, and non- decreasing, then there exists a unique solution for the 
fixed point equation ([7]) . 

Proof. The proof consists in showing that (jlSp has a unique solution. By the 
Implicit fmiction theorem, there exists a unique Xj(yj^i, j/j+i) such that, 

Xj{yj-i,yj+i) = + Xj{yj-i,yj+i),Xj{yj^i,y^+i)+yj+i), (16) 

and this function (yj-i, j/j+i) — s- Xj{yj-i,yj+i) is positive and continuous dif- 
ferentiable. Taking the partial derivative to on both sides of this identity, 
one gets that 



dxj{yj-i,yj+i) _ d(l>j{si,S2) 



dyj-i dsi 



^ [ ^^dxj{yj^i,yj+i) 



s=s{y) 



dyj-1 



d(t>j{si,S2) 



ds2 



^ dxj{yj-i,yj+i) 
dyj-i 
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with s{y) = {y.j-1 + x^- y^+i), y^+i) + yj+i). Hence, 



dy. 



'i-i 



;(gl,g2) 

dsi 



i>j{si,S2) d(t)j{si,S2) 



dsi 



ds2 



<0. 



s=s{y) 



A similar expression holds for dxj{yj^i, yj^i) / dyj-\-i < 0, and one can conclude 
that 



dxj{yj-i,yj+i) 



dyj~i 

' f d<j>j{si,S2) 



dxj{yj^i,y.j+i) 



dyj+i 

t>jisi,S2) 

ds2 



dsi 



■j{si,S2) 

ds2 



(17) 



8=s(y) 



If Xj(0,i/j+i) = 0, then by Relation (fT6|) . one gets that, for some constant aj, 
see Equation ([8]), 



= (l)j{Q,yj+i) = aj/^(3j{0,yj+i), 



which holds only if yj+i — oo. Hence, Xj{0,yj+i) > 0. Since Xj{yj^i,yj+i) is 
continuous and positive, and Xj(0, > 0, one obtains that there exists an 

M~ > such that yj-i + Xj{yj-i, yj+i) > MJ' for all yj+i > 0. Similarly, 
there exists an > such that Xj(j/j_i, j/j+i) + yj+i > . This gives the 
following upper bound, 

d(t)j{si,S2) 



dsi 



Oj d(3j{si,S2) 
2 



dsi 



< 



L 



2 (/3,(Mr,M+))3/2^ 



where we used that fUj is non-decreasing, Lipschitz continuous (with constant 
L) and differentiable. ^From Equation ([T7])) one now obtains that there exists a 
constant < C < 1 such that 



dxj{yj^i,yj+i) 



'3-1 



dxj{vj-i,yj+i) 



< C. 



Hence, the mapping T : M:[ with r(y) = (x^ (y^^i, j/j+i) for j = 1, . . . , J) 

is a contraction, and has a unique fixed point (j/p, i.e., Equation (|15p has a 
unique solution. 



Routes with one node or two consecutive nodes 

Consider now a ring with J nodes and K — 2J classes. Class j uses two nodes: 
nodes j and j + 1, j — 1, . . . , J. Class Oj uses one node: node j, j = 1, . . . , J. 
We assume that Aji^ = 1 if and only if class fc uses node j, and zero otherwise. 
We focus on functions that satisfy ([5]). Equation ([7]) is in this context 

Uj = 0oi(uj) + 4>3-i{u3-i + + <i>i{^3 + Uj+l), 1 < j < J- 
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Fig. 4. Routes with one node or two consecutive nodes 



For Hj = 4>j{uj + Uj+i) and j/oj = 4>aj{uj), j = 1 . . . , J, the above equation can 
be rewritten as foUows, for = 1, . . . , J, 



Proposition 7. // the functions Pk, l<k<J,Ol<k< OJ, are Lipschitz, 
continuously dijferentiable, non- decreasing, and satisfy (0), then there exists a 
unique solution for the fixed point equation ((T]). 

Proof. By the Imphcit function theorem, for each j, there exists a unique Xoj{t) 
satisfying the relation xoj{t) = (l)Qj{xoj(t)+t), and this function is non- increasing 
and continuous differentiable. One now has to solve the equation 

Vj = (t>]{y3-i + 2j/j + Vj+i + XQj{yj-i + yj) + xoj+i(yj + yj+i))- (19) 

^From the fact that —1 < x'^^^t) < 0, it can be easily checked that the right- 
hand side of Equation pop is non- increasing in yj. Hence, there exists a unique 
Xj{yj-i,yj+i) such that yj = Xj{yj^i,yj+i) satisfies Equation and this 

function yj+i) —f Xjlyj-iyyj+i) is positive and continuous differentiable 

(by the Implicit function theorem). In particular, Xj{yj^i, j/j+i) = 4>j{fj{y))-< for 
all j = 1, . . . , L, with 

fj{y) = Vj-i + 2a;j(yj_i, %+i) + yj+i 

+ ^OjiVj-i + Xj{yj-i,yj+i)) + xo^^+i{xj{yj^i,yj+i) + yj+i). 

^From this one can derive that 



Vj = ivj-i + 2% + Vj+i + yoj + yoj+i), 
yoj = <?!'oj(yoj +%■-! + j/j)- 



(18) 



dxj{yj_i,yj+i) 
dyj-i 



<l^'j{fAy)) [1 + Aj{yj-i + ^Ayj-i^yj+i))]/ 



1 - 4>'j{fjiy))('^ + 4i(%-i + xjiyj~iiyj+i)) 



Oj + l 



{xj{yj-i,yj+i) + yj+i)) < 0, 
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and a similar expression holds for dxj{yj-i,yj^i)/dyj+i < 0. As in the proof 
of Proposition [6l an upper bound on —<P'j{fj{y)) can be obtained. This implies 
that there exists a constant < C < 1 such that 



dxj{yj-i,yj+i) 



dyj- 



dxj{yj^i,yj+i) 



<C, J 



1, 



J. 



Hence, the mapping T : M![^ ^ R:[ with T{yi,...,yj) = {xj{yj-i,yj+i), j = 
1, . . . , J) is a contraction, and has a unique fixed point {y*j)- One concludes that 
there exists a unique solution y* and y'^^ — xoj{y*_i + yj), i = 1, of ^TE\i . 



Routes with two consecutive nodes and a complete route 

Consider a ring with J nodes and K = J + 1 classes. Class 1 < j < J uses two 
nodes: node j and j + 1 and class uses all nodes 1, . . . , J. 



Fig. 5. Routes with two consecutive nodes and a complete route 



We focus on fimctions (3k that satisfy (O. Equation (O is in this context 

Uj ~ <j)0 {Ul + ■ ■ ■ + Uj) + <j)j^l{Uj-l + Uj) + (j)j{Uj + Uj + l). 

For yj = 4>j{uj + and j/q = (poiui + U2 + ■ ■ ■ + uj), the above equation can 

be rewritten as follows 

Vj ^ (/^jiVj-i + '^yj +yj+i + 2yo), j = l,...,J, ^^^^ 
yo = 0o(>/2/o + 2 Y.j=i Vj)- 

Proposition 8. If the functions (3k, < k < J , are continuous, non- decreasing, 
satisfy and satisfy the assumptions of Proposition [H then there exists a 
unique solution for the fixed point equation ([7]). 

Proof. ^From Proposition [5] we have that Equation (|^D|) has at least one fixed 
point solution. Let x = (xj : j = 0, 1, . . . , J) and y — {yj : j = 0,1, J) both 
be fixed points. 
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If the relation yj < Xj holds for all j = 1, . . . , J, then 

(/»j(2?/j + + j/j+i + 2yo) = yj < Xj = (f)j{2xj + xj^i + Xj+i + 2xo). 

Since the function 0^ is non-increasing, one gets that 

2yj + y.j^i + yj+i + 2ya > 2xj + xj^i + Xj+i + 2xo. 

Summing over all = 1, . . . , L, we obtain 

4(yi + • ■ • + yj) + 2Jyo > 4(xi + • ■ • + xj) + 2Jxo, 

which implies that xq < j/o- However, 

yo = 00 (2(yi H \-yj) + Jyo) < (2(a;i H h a;,/) + Jxq) = xq, 

hence, we obtain a contradiction. 

We can conclude that there is an m G {1, . . . , J} such that j/m < Xm and 
J/m+i > Xm+i- To show that x — y, one proceeds along similar lines as in the 
proof of Proposition El 
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